TRANS-CANADA SLIMEWAYS: SLIME MOULD 
IMITATES THE CANADIAN TRANSPORT NETWORK 



ANDREW ADAMATZKY AND SELIM G. AKL 

Abstract. Slime mould Physarum polycephalum builds up so- 
phisticated networks to transport nutrients between distant part 
of its extended body. The slime mould's protoplasmic network is 
optimised for maximum coverage of nutrients yet minimum energy 
spent on transportation of the intra-cellular material. In labo- 
ratory experiments with P. polycephalum we represent Canadian 
major urban areas with rolled oats and inoculated slime mould in 
the Toronto area. The Plasmodium spans the urban areas with 
its network of protoplasmic tubes. We uncover similarities and 
differences between the protoplasmic network and the Canadian 
national highway network, analyse the networks in terms of prox- 
imity graphs and evaluate slime mould's network response to con- 
tamination. 

Keywords: biological networks, vehicular transport networks, slime 
mould, bio-inspired computing 



1. Introduction 

The increase of long-distance travel and subsequent reconfiguration 
of vehicular and social networks [14j requires novel and unconventional 
approaches towards analysis of dynamical processes in complex trans- 
port networks [7j, routing and localisation of vehicular networks [21j, 
optimisation of interactions between different parts of a transport net- 
work during scheduling road expansion and maintenance [26j , and shap- 
ing of transport network structure [8] . In the present paper we attempt 
to build viable analogies between biological and human-made transport 
networks and project behavioural traits of biological networks onto ex- 
isting vehicular transport networks. 

While choosing a biological object we want it to be experimental 
laboratory friendly, easy to cultivate and handle, and convenient to 
analyse its behaviour. Ants would indeed be the first candidate, and 
a great deal of impressive results has been published on ant-colony 
inspired computing [Wj [24], however ant colonies require substantial 
laboratory resources, experience and time in handling them. Actually 
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very few, if any, papers were published on experimental laboratory 
implementation of ant-based optimisation, the prevalent majority of 
publications being theoretical. There is however an object which is 
extremely easy to cultivate and handle, and which exhibits remarkably 
good foraging behaviour and development of transport networks. This 
is the Plasmodium of Physarum polycephalum. 

Plasmodium is a vegetative stage of acellular slime mould P. poly- 
cephalum^ a syncytium, that is, a single cell with many nuclei, which 
feeds on microscopic particles f25]. The Plasmodium is a unique user- 
friendly biological substrate from which experimental prototypes of 
massive-parallel amorphous biological computers are designed [4J . Dur- 
ing its foraging behaviour the Plasmodium spans scattered sources of 
nutrients with a network of protoplasmic tubes. The protoplasmic 
network is optimised to cover all sources of food and to provide a ro- 
bust and speedy transportation of nutrients and metabolites in the 
Plasmodium body. The Plasmodium's foraging behaviour can be in- 
terpreted as computation. Data are represented by spatial configura- 
tions of attractants and repellents, and results of computation by struc- 
tures of a protoplasmic network formed by the Plasmodium on the data 
sets [IHIIIHIII]. The problems solved by Plasmodium of P. polycephalum 
include shortest path |l6l^, implementatiton of storage modification 
machines [1] , Voronoi diagram [23j , Delaunay triangulation [4J , logical 
computing [29] , and process algebra ^22j ; see overview in [4J . 

Previously [2] we have evaluated a road-modeling potential of P. 
polycephalum^ however no conclusive results were presented back in 
2007. A step forward, namely, biological-approximation, or evalua- 
tion, of human-made road networks was done in our previous papers 
on approximation of motorways/highways in the United Kingdom [3j, 
Mexico [5j and the Netherlands [6] by Plasmodium of P. polycephalum. 
For all three countries we found that, in principle, the network of pro- 
toplasmic tubes developed by Plasmodium matches, at least partly, 
network of human-made transport arteries. The shape of a county and 
the exact spatial distribution of urban areas (represented by source 
of nutrients) may play a key role in determining the exact structure 
of the Plasmodium network. Also we suspect that a degree of match- 
ing between Physarum networks and motorway networks is determined 
by original government designs of motorways in any particular coun- 
try. This is why it is so important to collect data on development 
of Plasmodium networks in all major countries, and then undertake a 
comparative analysis. 

What are unique properties of the Canadian transport system? The 
Canadian Highway System gives us a good example of a logically 
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designed transportation system whose key goal is to connect all the 
provinces together by highways. The highway network was built as a 
federal-provincial territorial cooperative effort with great effort taken 
in coordinating work on different parts. Another attractive property 
of the highway system is that it was designed to provide an access to 
remote areas where no spots of high population density exist. 

The paper is structured as follows. We present experimental tech- 
niques used in Sect. [2} Properties of protoplasmic networks built by 
P. polycephalum are discussed in Sect. [3| We compare slime mould 
networks with the Canadian highway network in Sect. [4] and slime 
mould and human-made networks with proximity graphs in Sect. [5) In 
Sect. |6] we show how the slime mould transport network restructures 
in response to a spreading contamination. 



2. Methods 

The Plasmodium of P. polycephalum is cultivated in plastic contain- 
ers, on paper kitchen towels sprinkled with still drinking water and fed 
with oat flakes (Asda's Smart Price Porridge Oats). For experiments 
we use 12 x 12 cm polyestyrene square Petri dishes. Agar plates, 2% 
agar gel (Select agar, Sigma Aldrich), are cut in a shape of Canada. We 
consider the elven most populated urban areas U of Canada (Fig. [T^) 
and five transport nodes: 



(1) Toronto area (includ- 
ing Hamilton, London, 
St. CatharinesNiagara, 
Windsor, Oshawa, Barrie, 
Guelph, and Kingston) 

(2) Montreal area (incuding 
OttawaCatineau, Quebec 
City, Sherbrooke, Trois- 
Rivieres) 

(3) Vancouver area (includ- 
ing Victoria, Abbotsford, 
Kelowna) 

(4) Calgary 



(5) 


Edmonton 


(6) 


Winnipeg 


(7) 
(8) 


Halifax-Moncton 


Saskatoon- Regina 


(9) 


St. John's 


(10) 


Sudbury 


(11) 


Thunder Bay 


(12) 


Inuvik 


(13) 


Wrigley 


(14) 


YeUowknife 


(15) 


Thompson 


(16) 


Radisson 



The last five entries from Inuvik to Radisson, are not highly popu- 
lated urban areas. They are transport nodes added for completeness, 
i.e. to present slime mould with the same number of principle transport 
nodes as the human-made highways system (Fig. flk). Some transport 
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(b) 

Figure 1. Experimental setup, (a) Urban areas and 
transport nodes to be represented by oat flakes, from [9]. 
(b) Snapshot of protoplasmic transport network devel- 
oped by P. polycephum] the snapshot is made on a high- 
way map. 
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nodes as Fort McMurray, La Ronge, Flin Flon, and so on, are not in- 
cluded in the list due to their proximity to already chosen major urban 
areas. 

To project regions of U onto agar gel we place oat flakes in the 
positions of these regions of U (Fig. [T]). In each experiments we tried 
to match the size and shape of areas in (Fig. [T]) by selecting rolled oat 
of corresponding size and shape. At the beginning of each experiment 
a piece of Plasmodium, usually already attached to an oat flake in the 
cultivation box, is placed in the Toronto area (region 1 in Fig.[T^). The 
Petri dishes with Plasmodium are kept in darkness, at temperature 
22-25 C*^, except for observation and image recording. Periodically 
(usually in 12 h or 24 h intervals) the dishes are scanned in Epson 
Perfection 4490. We undertook 23 experiments. 



3. Foraging on urban areas 

It usually takes the Plasmodium of P. polycephalum 2-5 days to span 
all urban areas. How fast the Plasmodium colonises the space depends 
on many unknown factors, including seasonal variations, Plasmodium's 
age, etc. 'Younger' plasmodia, which were just recently 'woken up' 
from the sclerotium phase do usually colonise the experimental arena 
quicker than old plasmodia, which were replanted several times in cul- 
ture boxes. Images of protoplasmic networks presented in the paper 
are taken when all oat flakes, representing U, were colonised by Plas- 
modium. Examples of the protoplasmic networks are shown in Fig. [2j 

As every living creature does, the Plasmodium of P. polycephalum 
rarely repeats its foraging pattern, and almost never builds exactly 
the same protoplasmic network twice. To generalise our experimen- 
tal results we constructed a Physarum graph with weighted-edges. A 
Physarum graph is a tuple P = (U,E,t(;), where U is a set of urban 
areas, E is a set edges, and w : E ^ [0^1] associates each edge of E 
with a probability (or weights). For every two regions a and b from 
U there is an edge connecting a and 6 if a Plasmodium's protoplasmic 
link is recorded at least in one of k experiments, and the edge (a, b) 
has a probability calculated as a ratio of experiments where protoplas- 
mic link (a, b) occurred in the total number of experiments k = 23. 
For example, if we observed a protoplasmic tube connecting areas a 
and 6 in 5 experiments, the weight of edge (a, 6) will be w{a^b) = ^. 
We do not take into account the exact conflguration of the protoplas- 
mic tubes but merely their existence. Further we will be dealing with 
threshold Physarum graphs P(^) = (U, T(E), 0). The threshold 
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Figure 2. Examples of protoplasmic networks devel- 
oped by P. polycephalum on major urban areas and 
transport nodes U. Grey-scale images. 
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Figure 3. Configurations Physarum-gidiph P{0) for 
various cutoff values of 9. Thickness of each edge is pro- 
portional to the edge's weight. 
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Physarum graph is obtained from Physarum graph by the transforma- 
tion: T(E) = {e G E : w{e) > 9}. That is aU edges with weights 
less than or equal to 9 are removed. Examples of threshold Physarum 
graphs for various values of 9 are shown in Fig. [3| 

A 'raw' Physarum graph P(0) is a non-planar graph due to the pres- 
ence of protoplasmic tube connecting Saskatoon-Regina (8) with Yel- 
lowknife (14) (Fig. |3^). It also exhibits two cross-Canada transport 
links Inuvik (12) to Radisson (16) and Sudbury (10) to Yellowknife 

(14) . Nevertheless, these are links that might be considered as sense- 
less from a geographical point of view because they are crossing massive 
of mountains and forests. 

All three links disappear when we increase ^ to 8 (Fig. ^p). Four 
more links become trimmed off when 9 = 8: Sudbury (10) — Thompson 

(15) , Montreal area (2) — Sudbury (10), Sudbury (10) — Radisson 

(16) , Montreal area (2) — St. John's (9). 

P(^) is the last connected graph in a series of threshold Physarum 
graphs, P(^), < ^ < 23. Urban areas St. Johns's (9) and Radisson 
(16) become isolated in P(^) due to the disappearance of transport 
links from the Montreal area (2) to Radisson (16) and from Radisson 
to St. John's (9) (Fig.|3|3). 

The Physarum graph splits into a tree and two isolated nodes for 9 = 
II (Fig. pp). Five long distance routes from urban areas to transport 
nodes are removed: Vanvouver area (3) - Inuvik (12), Vancouver area 
— Wrigley (13), Edmonton (5) — Wrigley, Saskatoon-Regina (8) — 
Thompson (15), Winnipeg (6) — Thompson; and, two short-distance 
routes between major urban areas: Vancouver area (3) — Edmonton 
(5) and Calgary (4) — Saskatoon-Regina (8). 

The graph P(||) consists of several disconnected components: two 
isolated nodes, St. John's and Sudbury, two-node segment Inuvik — 
Wrigley, and a tree spanning the rest of the urban ares (Fig. [3^). A 
further increase of ^ to || (Fig. |3]f) leads to the formation of 

• five isolated nodes: Vancouver area, St. John's, Yellowknife, 
Thompson, Winnipeg; 

• a segment: Inuvik — Wrigley; 

• a chain: Calgary — Edmonton — Saskatoon-Regina — Win- 
nipeg — Thunder Bay; 

• a chain: Sudbury — Toronto area — Montreal area — Halifax- 
Moncton. 

Only segment-routes Calgary to Edmonton and Winnipeg to Thun- 
der Bay are represented by protoplasmic tubes in almost all experi- 
ments (Fig. |3^). 
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(e) P(i)nH 

Figure 4. Physamm vs highway network, (a) Highway 
network in Canada [9j. (b) Highway graph H. (cde) In- 
tersection of threshold Physarum graph with highway 
graph for ^ = 0, ^, 
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4. Physarum network vs highway network 



We construct the highway graph H as foUows. Let U be a set of 
urban regions, for any two regions a and b from U, the nodes a and 
b are connected by an edge (a, b) if there is a motorway starting in 
the vicinity of a and passing in the vicinity of b and not passing in the 
vicinity of any other urban area c G U. Highway graph H shown in 
Fig. [4b is extracted from a scheme of the Canadian transport network 



in the vicinity of a and passing in the vicinity of b and not passing 
in the vicinity 

Finding 1. Physarum almost approximates the Canadian highway net- 
work. 

'Raw' Physarum graph P(0) approximates 21 of 22 edges of highway 
graph H. Only one edge, Vancouver area to Calgary, of the highway 
graph H is not represented by protoplasmic tubes in any of the 23 
experiments undertaken (Fig. The 'raw' physarum graph gives us 
a rather relaxed approximation because it includes even links which 
occurred just once in a set of experiments. Let us look at P(^) which 
represents links which occurred in over 35% of experiments. Physarum 
graph P(^) approximates 18 of 22 edges of H (Fig. [ijl). The only 
edges of H not represented in P(^) are Vancouver area to Yellowknife, 
Sudbury to Radisson, Sudbury to Montreal area, and St. John's to 
Montreal area. 

Finding 2. A core component of the Physarum transport network and 
the Canadian highway network consists of a chain passing along the 
south border from Halifax-Monctron area to Edmonton^ and a fork 
attached to Edmonton; the south branch of the fork is Edmonton — 
Calgary — Vancouver area and the north branch is Edmonton - Yel- 
lowknife - Wrigley. 

The component above is the only connected component in the inter- 



A planar graph consists of nodes which are points of the Euclidean 
plane and edges which are straight segments connecting the points. A 
planar proximity graph is a planar graph where two points are con- 
nected by an edge if they are close in some sense. A pair of points is 
assigned a certain neighbourhood, and points of the pair are connected 




section of P(||) with H (Fig. |4^). 



5. Proximity graphs 
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(a) GG (b) RNG = MST 




(c) GGpH (d) RNGpH 



Figure 5. Proximity graphs and their intersection with 
highway graph H. (a) Gabriel graph, (b) Relative neigh- 
bourhood graph and minimun spanning tree, (c) Inter- 
section of Gabriel graph with highway graph, (d) In- 
tersection of relative neighbourhood graph with highway 
graph. 

by an edge if their neighbourhood is empty. Here we consider the most 
common proximity graph as follows. 

• GG: Points a and b are connected by an edge in the Gabriel 
Graph GG if disc with diameter dist{a^ h) centered in middle 
of the segment ah is empty ^iTl [15] (Fig. [5^). 

• RNG: Points a and h are connected by an edge in the Relative 
Neighbourhood Graph RNG if no other point c is closer to a 
and h than dist{a, h) [28j (Fig. |5|3). 

• MST: The Euclidean minimum spanning tree (MST) [20j is 
a connected acyclic graph which has minimum possible sum of 
edges' lengths (Fig. [5|3). 
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In general, the graphs relate as MST C RNG C GG [28l [TSj [T2] . 

Finding 3. For a given configuration of nodes ofU RNG = MST. 

The finding implies that the configuration of urban areas of U is 
'spanning friendly'. 

Finding 4. Let T^ be a minimum spanning tree rooted in node i G U 
then Ti = Tj for any i^j G U. 

We demonstrated this by direct computation of all possible spanning 
trees on U. 

Finding 5. St. John^s and Inuvik urban areas are isolated in GG f] RNG 
and GGHMST. 

Finding 6. NIST /{(Inuvik — Wrigley), (Halifax-Moncton — St. John^s)} 
C H 

This means that Canadian highway network is almost optimal (Fig.[5]le). 

Intersections of Physarum graphs for principle values of threshold 9 
with the Gabriel graph and the minimum spanning tree are shown in 
Fig.|6| 

Finding 7. GG /{(Inuvik — Wrigley), (Thunder Bay — Radisson)} 
C H 

Finding 8. MST C P(0) 

'Raw' Physarum includes an 'ideal' acyclic spanning network MST. 
This somehow characterises a good quality of a slime mould approx- 
imation of a transport network. The minimum spanning tree is not 
included in the high-threshold Physarum graph P(||). However, there 
is a 'strong' component of MST which is included in the high-threshold 
Physarum graph (Fig. [6]). The strong component is a tree rooted in the 
Toronto area. The tree's stem is Toronto — Winnipeg — Saskatoon- 
Regina — Edmonton. It has three branches. A small branch Winnipeg 
— Thompson, and two end branches Edmonton — Calgary — Vancou- 
ver area and Edmonton — Yellowknife — Wrigley — Inuvik. 

6. Response to contamination 

To imitate propagating contamination we place a crystal of sea salt 
(SAXA Coarse Sea Salt, a crystal weight around 20 mg), in the place of 
Bruce Nuclear Power station. Inorganic salts are chemo-repellents for 
P. polycephalum therefore sodium chloride diffusing in agar gel causes 
Plasmodium to retreat from a contaminated zone. We studied Plas- 
modium's response circa 24 h after initiation of contamination. During 
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(e)P(i)nGG (f)P(i)nMST 

Figure 6. Intersection of (ace) Gabriel graph GG and 
(bdf) minimum spanning tree MST with Physarum 
graphs P(^) for (ab) ^ = 0, (cd) 9^^, (ef) 9 = ^. 



24 h a contamination zone expands as far as Winnipeg an Thompson 
in the west and St. John's in the east. In some cases contamination 
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(e) (f) 

Figure 7. Physarum response to contamination, 
(ab) Migration outside Canada, (c-f) Compensating ac- 
tivation of transport network unaffected by contamina- 
tion, (ace) Original images, (bdf) Binarized images: 
only pixels from (ace) which red and green components 
exceed 100 and blue component is less than 100 are 
drawn as black pixels in (bdf), otherwise white. 
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(c) (d) 

Figure 8. Physarum sprouting in response to contami- 
nation, (ab) Sprouting from edges, (cd) Sprouting from 
nodes, (ac) Original colour images, (bd) Binarized im- 
ages. 



spreads till Saskatoon-Moncton. In a few experiments Plasmodium 
colony occupying St. Jonn's remains unaffected. 

Finding 9. In response to contamination propagating from Bruce Nu- 
clear Power station^ the Plasmodium of P. polycephalum takes one or 
more of the following actions: migrates outside Canada^ enhances the 
transport network outside the contaminated zone, sproutes indiscrimi- 
nately from urban areas and transport links. 

The Plasmodium's reactions are illustrated in Figs. [7| and [8| Four 
types of responses are observed in laboratory experiments. 
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• Plasmodium migrates outside Canada (Fig.[7^b). Typical waves 
of migration are from Nunavut towards Baffin Bay and Green- 
land and from British Columbia towards Washington and Ore- 
gon in USA. Due to the growth substrate being an agar plate 
cut in the shape of Canada, the Plasmodium ends up on the bare 
plastic bottom of a Petri dish, therefore it does not migrate far 
away from Canada. 

• Plasmodium enhances it foraging and colonisation activity in 
the parts unaffected by contamination (Fig. [Tj^-f), mainly in 
Alberta, British Columbia, Northern Territories, Yukon and 
Nunavut. Protoplasmic tubes themselves are often increased in 
size and intensity of their colours, which reflects increased prop- 
agation of cytoplasm inside the tubes. For example, in Fig. 
Plasmodium clearly shows hyper-activity in the Nunavut area, 
with the whole territory covered by spreading Plasmodium. Fig. 
illustrates hyper-activation of transport routes, particularly links 
Inuvik — Wrigley, Yelloknife — Wrigley, Edmonton — Wrigley, 
Edmonton — Vancouver area. Saskatoon- Regina — Edmonton, 
and Winnipeg — Saskatoon. 

• Plasmodium expands outside oat flakes and also produces pro- 
cesses protruding from the protoplasmic tubes (Fig.jS^-d). This 
a common reaction of Plasmodium in response to mechanical 
damage, e.g. cutting of protoplasmic tubes [4J. 

7. Discussion 

To imitate transport networks in Canada we represented major ur- 
ban areas and transport nodes with oat flakes, inoculated Plasmodium 
of Physarum polycephalum^ allowed the Plasmodium to span all oat 
flakes with its network of protoplasmic tubes and analysed the struc- 
ture of the protoplasmic network. We found that in over 75% of ex- 
periments P. polycephalum an acyclic transport network consisting of 
a chain spanning urban areas along the south boundary of Canada, 
from Halifax-Moncton to Vancouver area, with branches Edmonton to 
Yellowknife to Wrigley to Inuvik and Yellowknife to Thompson. In 
all experiments slime mould approximates all but Vancouver to Cal- 
gary links of Canadian highways networks. Both slime mould and 
Canadian highway networks have a strong spanning tree component 
and thus can be thought of as optimal transport networks. In labora- 
tory experiments with slime mould we also detailed possible scenarios 
of transport network restructuring in a response to a spreading con- 
tamination. We believe our results make a substantial contribution 
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towards nature-inspired analyses and design of human-made transport 
networks. Further experiments are necessary to determine how nat- 
ural and geographical conditions, especially terrain, affect the exact 
topology of developing transport networks. 
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